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1 Introduction 



The purpose of this paper is to extend classical stochastic calculus for multi- 
dimensional Brownian motion to the setting of nonlinear G-expectation. We 
first recall the general framework of nonlinear expectation studied in [44j and 
[43] , where the usual linearity is replaced by positive homogeneity and subaddi- 
tivity. Such a sublinear expectation functional enables us to construct a Banach 
space, similar to an L 1 -space, starting from a functional lattice of DanielPs type. 

Then we proceed to construct a sublinear expectation on the space of contin- 
uous paths from M.+ to R d , starting from 0, which will be an analogue of Wiener's 
law. The operation mainly consists in replacing the Brownian semigroup by a 
nonlinear semigroup coming from the solution of a nonlinear parabolic partial 
differential equation |T]) where appears a mapping G acting on Hessian matri- 
ces. Indeed, the Markov property permits to define in the same way nonlinear 
conditional expectations with respect to the past. Then we presents some rules 
and examples of computations under the newly constructed G-Brownian (mo- 
tion) expectation. The fact that the underlying marginal nonlinear expectations 
are G-normal distributions derived from the nonlinear heat equation |T]) is very 
helpful to estimate natural functionals. As result, our G-Brownian motion also 
has independent increments with identical G-normal distributions. 

G-Brownian motion has a very rich and interesting new structure which non 
trivially generalizes the classical one. We thus can establish the related stochas- 
tic calculus, especially G-Ito's integrals (see 1942]) and the related quadratic 
variation process (B) . A very interesting new phenomenon of our G-Brownian 
motion is that its quadratic process (B) also has independent increments which 
are identically distributed. The corresponding G-Ito's formula is obtained. We 
then introduce the notion of G-martingales and the related Jensen inequality 
for a new type of "G-convex" functions. We have also established the exis- 
tence and uniqueness of solution to stochastic differential equation under our 
stochastic calculus by the same Picard iterations as in the classical situation. 
Books on stochastic calculus e.g., [10], [23], [25], [27], [31], [36], [50], [5T], [SS] are 
recommended for understanding the present results and some further possible 
developments of this new stochastic calculus. 

As indicated in Remark[5] the nonlinear expectations discussed in this paper 
can be regarded as coherent risk measures. This with the related conditional 
expectations E[-|7i t ]t>o makes a dynamic risk measure: G-risk measure. 

The other motivation of our G-expectation is the notion of (nonlinear) g— 
expectations introduced in [35], [JD]. Here g is the generating function of a 
backward stochastic differential equation (BSDE) on a given probability space 
(Q, J 7 , P). The natural definition of the conditional ^-expectations with respect 
to the past induces rich properties of nonlinear ^-martingale theory (see among 
others, [3J, [5], @, [7], [11], [12], [8], [9], [21], [29], [41], [45], [46], [48]). Recently 
^-expectations are also studied as dynamic risk measures: (/-risk measure (cf. 
[52] . [4], [17]). Fully nonlinear super-hedging is also a possible application (cf. 
[3"5] and [53] where new BSDE approach was introduced). 

The notion of g-expectation is defined on a given probability space. In [44j 



2 



(see also [13]), we have constructed a kind of filtration-consistent nonlinear 
expectations through the so-called nonlinear Markov chain. As compared with 
the framework of (/-expectations, the theory of G-expectation is intrinsic, a 
meaning similar to "intrinsic geometry" in the sense that it is not based on a 
given (linear) probability space. Since the classical Brownian expectation as 
well as many other linear and nonlinear expectations are dominated by our G— 
Expectation (see Remark[25l Example l40l and [44] ) and thus can be considered as 
continuous functionals, our theory also provides a flexible theoretical framework. 

1-dimensional G-Brownian motion was studied in [J7] . Unlike the classical 
situation, in general, we cannot find a system of coordinates under which the 
corresponding components B % , i = 1, • • • , d, are mutually independent from each 
others. The mutual quadratic variations (B l , B- 5 ) will play essentially important 
rule. 

During the reversion process of this paper, the author has found a very in- 
teresting paper [18] by Denis and Martini on super-pricing of contingent claims 
under model uncertainty of volatility. They have introduced a norm on the 
space of continuous paths = G([0,T]) which corresponds to our L^-norm 
and developed a stochastic integral. There is no notions of nonlinear expecta- 
tion such as G-expectation, conditional G-expectation, the related G-normal 
distribution and the notion of independence in their paper. But on the other 
hand, powerful tools in capacity theory enables them to obtain pathwise results 
of random variables and stochastic processes through the language of "quasi- 
surely", (see Feyel and de La Pradelle [21]) in the place of "almost surely" in 
classical probability theory. Their method provides a way to proceed a pathwise 
analysis for our G-Brownian motion and the related stochastic calculus under 
G-expectation, see our forthcoming paper joint with Denis. 

This paper is organized as follows: in Section 2, we recall the framework of 
nonlinear expectation established in (44] and adapt it to our objective. In sec- 
tion 3 we introduce e?-dimensional G-normal distribution and discuss its main 
properties. In Section 4 we introduce d-dimensional G-Brownian motion, the 
corresponding G-expectation and their main properties. We then can estab- 
lish stochastic integral with respect to G-Brownian motion of Ito's type, the 
related quadratic variation processes and then G-Ito's formula in Section 5, G- 
martingale and the Jensen's inequality for G-convex functions in Section 6 and 
the existence and uniqueness theorem of SDE driven by G-Brownian motion in 
Section 7. 

The whole results of this paper are based on the very basic knowledge of 
Banach space and the parabolic partial differential equation (TT]). When this 
G-heat equation (fTJ) is linear, our G-Brownian motion becomes the classical 
Brownian motion. This paper still provides an analytic shortcut to reach the 
sophistic Ito's calculus. 
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2 Nonlinear expectation: a general framework 



We briefly recall the notion of nonlinear expectations introduced in [44] . Follow- 
ing Daniell's famous integration (cf. Daniell 1918 [14] . see also [54]), we begin 
with a vector lattice. Let 17 be a given set and let TL be a vector lattice of real 
functions defined on Q containing 1 , namely, TL is a linear space such that 1 € TL 
and that X 6 TL implies \X\ £ TL. TL is a space of random variables. We assume 
the functions on TL are all bounded. 

Definition 1 A nonlinear expectation E is a functional TL i— > M satisfying 
the following properties 

(a) Monotonicity: if X,Y e TL and X >Y then E[X] >E[Y]. 

(b) Preserving of constants: E[c] = c. 

in this paper we are interested in the sublinear expectations which satisfy 

(c) Sub-additivity (or self dominated property): 

E[X]-E\Y] <E[X-Y], VX,YeTL. 

(d) Positive homogeneity: E[XX] = XE[X], VA > 0, X G TL. 

(e) Constant translatability: E[X + c] = E[X] + c. 

Remark 2 It is clear that (d) + (e) implies (b). We recall the notion of the 
above sublinear expectations was systematically introduced by Artzner, Delbaen, 
Eber and Heath JBj, in the case where £1 is a finite set, and by Delbaen 
in general situation with the notation of risk measure: p(X) = E[— X] . See also 
in Huber \2^j for even early study of this notion E ( called upper expectation E* 
in Ch.10 ofJEfl). 

We follow [33] to introduce a Banach space via TL and E. We denote ||X|| := 
E[|X|], X £ TL. TL forms a normed space (TL, ||-||) under ||-|| in the following 
sense. For each X, Y e TL such that \\X - Y\\ = 0, we set X = Y. This is 
equivalent to say that the linear subspace 

Ho :={XeTL, \\X\\ =0} 

is the null space, or in other words, we only consider the elements in the quotient 
space TL/TLq. Under such arrangement (TL, ||-||) is a normed space. We denote by 
([TL], ||-||), or simply [TL], the completion of (TL, ||-||). (TL, ||-||) is a dense subspace 
of the Banach space ([TL], ||-||) (see e.g., Yosida [56] Sec. 1-10). 
For any X € TL, the mappings 

X+(uj) = max{Z(cj),0} :H — ► TL, 
X~(u) = max{-X(o;),0} : TL i — > TL 
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satisfy 



\X+ -Y+\ < \X-Y\, 
X~ - Y~ <(Y- X) + < \X - Y\. 

Thus they are both contract mappings under ||-|| and can be continuously ex- 
tended to the Banach space [Tt] . 

We define the partial order ">" in this Banach space. 

Definition 3 An element X in ([H], ||-||) is said to be nonnegative, or X > 0, 
0<X,ifX = X+. We also denote by X >Y , orY <X. if X - Y > 0. 

It is easy to check that if X > Y and Y > X, then X = Y in ([H], ||-||). The 
nonlinear expectation E[-] can be continuously extended to ([H], ||-||) on which 
(a)-(e) still hold. 

3 G— normal distributions 

For a given positive integer n, we will denote by (x, y) the scalar product of x, 
y € R n and by \x\ — (x, x) 1 ! 2 the Euclidean norm of x. We denote by lip{R n ) the 
space of all bounded and Lipschitz real functions on R™. We introduce the notion 
of nonlinear distribution- G-normal distribution. A G— normal distribution 
is a nonlinear expectation defined on lip(M. d ) (here R d is considered as ft and 
lip(R d ) as H): 

P 1 G (0) = u(l, 0) : e lip(M d ) i-> R 

where u = u(t, x) is a bounded continuous function on [0, oo) x R d which is 
the viscosity solution of the following nonlinear parabolic partial differential 
equation (PDE) 

du 

— -G{D 2 u) = 1 u(0,x) = <f>(x), (t,x) e [0,oo) x R d , (1) 
here D 2 u is the Hessian matrix of u, i.e., D 2 u = {d* ixJ u)f ,- =1 and 

G(A) = G r (A) = \ suptr[ 77 T ^], A = (^)?j=i e S d . (2) 
z ~tev 

§d denotes the space of d x d symmetric matrices. T is a given non empty, 
bounded and closed subset of R dxd , the space of all d x d matrices. 

Remark 4 The nonlinear heat equation {!]] is a special kind of Hamilton- 
Jacobi- Bellman equation. The existence and uniqueness of (QJ) in the sense 
of viscosity solution can be found in, for example, J13\/ . J2zj, f38\j . J55f . and i32\/ 
for C 1 ' 2 -solution if"fj T > cr^In, for each 7 6 T , for a given constant cto > (see 
also in f36f for elliptic cases). It is a known result that u(t, •) £ lip(R d ) (see e.g. 
J55]j Ch.4, prop. 3.1. or ' L 38l Lemma 3.1. for the Lipschitz continuity ofu(t,-), 
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or Lemma 5.5 and Proposition 5.6 in '4-31 for a more general conclusion). The 
boundedness is simply from the comparison theorem (or maximum principle) of 
this PDE. It is also easy to check that, for a given ip £ Hp(R d xR d ), PfOO, •)) 
is still a bounded and Lipschitz function in x. 

In the case where T is a singleton {70}, the above PDE becomes a stan- 
dard linear heat equation and thus, for G° = G{ 70 }, the corresponding G°- 
distribution is just the d-dimensional classical normal distribution A/"(0, 7o7(f )■ 
In a typical case where 70 = Id £ I\ we have 

In the case where 70 £ T, from comparison theorem of PDE, 

FfW^fW.V^MK"). (3) 

More generally, for each subset V C T, the corresponding P Gr ' -distribution is 
dominated by P G in the following sense: 

Pf r ' (0) - Pf r ' (VO < P?{<p - VO, S lip(R d ). 

Remark 5 In '^71 we have discussed 1-dimensional case, which corresponds 
d = 1 and T = [a, 1] C R, where £ [0, 1] is a given constant. In this case the 
nonlinear heat equation (QJ) becomes 

^- ^[{dl x u) + -a 2 {dl x u)-]=Q, u(Q,x) = <f>(x), (M)6[0 lM )xI. 

In multi- dimensional case we also have the following typical nonlinear heat equa- 
tion: 

i—l 

where o~i £ [0, 1] are given constants. In this case we have 

T = {diag[ji, ■ ■ ■ ,7d],7» S [0*, 1], i = I,--- ,d}. 

The corresponding normal distribution with mean at x £ M. d and square 
variation t > is P^((f)(x + \ft x •)). Just like the classical situation of a normal 
distribution,, we have 

Lemma 6 For each <\> £ lip(M. d ), the function 

u(t,x) = P^{(j)(x + Vtx-)), (t, x) £ [0, 00) x R d (4) 

is the solution of the nonlinear heat equation {Ip with the initial condition 
u(0, •) = #•)• 
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Proof. Let u G C([0,oo) x R d ) be the viscosity solution of (P) with u(0, •) = 
</>(■) G lip(R d ). For a fixed (t,x) € (0, oo) x R d , we denote u(t,x) = u(t x 
i,xVt + x). Then u is the viscosity solution of JT]) with the initial condition 
u(0, a;) = 4>(xVt + x). Indeed, let ip be a C 1,2 function on (0, oo) x M d such that 
ip > u (resp. ip < u) and iP(t 7 £) — u(t, £) for a fixed (r, £) € (0, oo) x R d . We 
have ?/>(f, > u(t,x), for all (i,x) and 

^ x _ ^ ^ ^ ^ _ ^ _^ 

Since it is the viscosity solution of |T]), at the point (t, x) = (ri^Vt + x), we 
have 

But G is a positive homogenous function, i.e., G(\A) — XG(A), when A > 0, we 
thus derive 

^t^l-G(D 2 ij(t,x))\ (t ^ ={T ^<0 (resp. > 0). 

This implies that u is the viscosity subsolution (resp. supersolution) of (Q]). 
According to the definition of P (•) we obtain ■ 

Definition 7 FFe denote 

Pf((p)(x) = Pf(<p(x + Vix ■)) =u(t,x), (t, x) G [0, oo) x M d . (5) 

From the above lemma, for each <p G lip(M. d ), we have the following nonlinear 
version of chain rule: 

P t G (P?(<P))(x) = Pf +s W(x), Me [0,oo), ^ G K d . (6) 

This chain rule was initialled by Nisio [34] and [35] in terms of "envelope of 
Markovian semi-groups" . See also [44] . 

Lemma 8 The solution of {!]) with initial condition u(0,x) = 0((a, x)), for a 
given (p G lip(R), has the form u(t,x) = u(t,x), x = (a, x), where u is the 
solution of 

Ou 

-57 - G^d^u) = 0, u(0,x)=<p(x), (t,x)€ [O.ooJxR, (7) 

w/iere 

G a (/3) = i max tr [7 7 T aa T /3] , /3 G M. 
2 7er 

TTie above PDE can be written 
du i 

07 - 2 [ CT aa^ (9sxU) + + (T_ aa T (c^u) - ] = 0, u(0, S) = 0(aj). (8) 



7 



where we denote aa T = [a 1 a^] d G §d o,nd 

a A = sup tr [jj T A] = 2G(A), AeS d . (9) 

Here Ed is the space of d x d symmetric matrices. 

Remark 9 It is clear that the functional 

Pf«((p) = tZ(l,0) : € lip(R) h-> R 

constitutes a special 1-dimensional nonlinear normal distribution, called G a - 
normal distribution. 

Proof. It is clear that the PDE has a unique viscosity solution. We then 
can set u(t, x) = u(t, (a, a:)) and check that u is the viscosity solution of (JTJ. (JSJ) 
is then easy to check. ■ 

Example 10 In the above lemma, if <f> is convex, and o~ aa T > 0, then 

P t G (cj>(( a ,.))(x) = / 0(y)exp(-^ L)dy. 

If (j) is concave and cr_ aa T < 0, then 

P t G (^((a,.))(x) = / ^Jexpf-^ Udy. 

Proposition 11 We have 

(i) For each t > 0, the G -normal distribution P t G is a nonlinear expectation on 
the lattice lip(M ), with = M. d , satisfying (a)-(e) of definition^ The corre- 
sponding completion space [H] = [lip(Bi d )]t under the norm \\4>\\ t := P t G (|0|)(O) 
contains 4>(x) = a;™ 1 x • ■ ■ x x n d d , n,- = 1, 2, • • • , i = 1, • • • , d, x = (xi, ■ ■ ■ , Xd) T 
as well as x\ n x • • ■ x x n d d x tp(x), ip € lip(R d ) as its special elements. Relation 
(5p still holds. We also have the following properties 

(ii) We have, for each a = (ax, ■ • • , a<j) T € M. d and A £ E>d 

Pf (( a ' x )x£R d ) = 0, 

p t C ((( a > a; ) 2 )xeR«) = t ■ <W, P G ((-(a,x) 2 ) xeRd ) = t ■ o_ aa T, 
P t G (((a,x) 4 ) xmd ) = 6(a aa T)t 2 , P t G ((-(a,a;) 4 ) x6Ii<i ) = -6(a_ aaT ) 2 t 2 , 
P t G (((Ax, x)) xmd ) =t-a A = 2G(A)t. 

Proof, (ii) By Lemma [8j we have the explicit solutions of the nonlinear PDE 
([l} with the following different initial condition it(0, x) = 4>(x): 

4>{x) = (a, a;) ==> u(t, x) = (a, ar), 

4>(x) = (a, a;) 4 =>• u(t, x) = (a, a;) 4 + 6 (a, x) 2 a a3 rt + 6o- 2 aaT t 2 , 
<f>(x) — —(a, a:) 4 ==> u(t, x) = — (a, a;) 4 + 6 (a, x) 2 a_ aa Tt — 6|u_ aa r \ 2 t 2 . 
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Similarly, we can check that cj)(x) — (Ax, x) =^> u(t, x) = (Ax, x) + aAt. This 
implies, by setting A = aa T and A = — aa T , 

4*(x) = (a, x) 2 => u(t, x) = (a, x) 2 + (T &a rt, 
<p(x) = —(a, x) 2 =^> u(t, x) = —(a, x) 2 + a_ aa Tt. 

More generally, for 4>(x) = (a,x) 2n , we have 

u(t,x) = -j=L= r y 2n eM- { i^)dy. 
By this we can prove (i). ■ 

4 G— Brownian motions under G— expectations 

In the rest of this paper, we set ft = C d (R + ) the space of all R d -valued contin- 
uous paths (a>t)t£M+, with ujq — 0, equipped with the distance 

oo 

p(u\u 2 ) :=^2-*[(max \cu} ~ uj 2 \) A 1]. 

i—l ' 

17 is the classical canonical space and u> = (u t ) t >o is the corresponding canonical 
process. It is well-known that in this canonical space there exists a Wiener mea- 
sure (fi, T, P) under which the canonical process B t (uj) — u t is a d-dimensional 
Brownian motion. 

For each fixed T > 0, we consider the following space of random variables: 

L%(H T ) :={X( W ) = 9 iK,-, Wi J,Vm>l, h,---, t m G [0, T], e h P (R dxm )}- 

It is clear that {L® p (H.t)}t>o constitute a family of sub-lattices such that L^ p (Ti. t ) C 
L® p (TLt), for t < T < oo. L® p (TL t ) representing the past history of cu at the time 
t. It's completion will play the same role of Brownian filtration Tf as in classical 
stochastic analysis. We also denote 

oo 

L\ p (U) := \jL%(H n ). 

n=l 

Remark 12 It is clear that lip(R dxm ) and then L° p (H T ), L° lp (H) are vector 
lattices. Moreover, since 4>,ip G lip(R dxm ) implies <j) ■ ip G lip(R dxm ) thus X, 
Y G L%(H T ) implies X ■ Y G L%(H T ); X, Y G L%(H) implies X ■ Y G L a lp (H). 

We will consider the canonical space and set B t (uj) — u>t, t G [0,oo), for 
u G fi. 

Definition 13 The canonical process B is called a (d- dimensional) G-Brownian 
motion under a nonlinear expectation E defined on L\ '(H) if 
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(i) For each s,t > and ip E lip(M. d ), B t and B t + S — B s are identically dis- 
tributed: 

E[^(B t+s - B s )\ = E[ib(B t )] = if ty). 

(ii) For each m = 1, 2, • • • , < t\ < • • • < t m < oo, the increment B tm — B tm _ 1 
is "backwardly" independent from B tl ,- • ■ 1 B tm _ 1 in the following sense: for 
each cj> € lip(R dxm ), 

E&iBt,,--- ,B tm _ 1 ,B tm )]=E[MB tl ,--- ,B tm ^)], 

whereMx 1 ,--- ,x m ~ 1 ) = E[^(x 1 , • • • , x™- 1 , B^-B^^+x" 1 ' 1 )}, x 1 ,■ ■ ■ ,x m - 1 € 
R d . 

The related conditional expectation of 4>(B tl , • • • , B tm ) under H. tk is defined by 
E[0(B tl ,--- ,B tk ,--- ,B tm )\H tk } = <f> m - k (B tl ,'-- ,B tk ), (10) 

where 

fim-kix 1 , ■ ■ ■ ,x k )=E[<f>(x 1 ,--- ,x k ,B tk+1 - B tk +x k ,--- ,B tm - B tk +x k )}. 

It is proved in [33] that E[-] consistently defines a nonlinear expectation on 
the vector lattice L® p (Ht) as well as on L® p (H) satisfying (a)-(e) in Definition 
[TJ It follows that E[|X|], X e L° p (H T ) (resp. L° lp {H)) forms a norm and thus 
L° ip (H.T) (resp. L^ p (TL)) can be extended, under this norm, to a Banach space. 
We denote this space by L g (Ht) (resp. L G (Tt)). For each < t < T < oo, 
we have L l G (H t ) C L g (Ht) C L g (H). In L G (H T ) (resp. L G (H T )), E[-] still 
satisfies (a)-(e) in Definition [TJ 

Remark 14 It is suggestive to denote L® p (H. t ) by 7i® and L G (Tlt) by Ht, 
L G (7i) by Tt and thus consider the conditional expectation E[-|7i t ] as a pro- 
jective mapping from 7i to Tit- The notation L G (Tlt) is due to the similarity of 
L (fl,!Ft,P) in classical stochastic analysis. 

Definition 15 The expectation E[-] : L G (7i) i— > K introduced through above 
procedure is called G— expectation, or G~Brownian expectation. The corre- 
sponding canonical process B is said to be a G-Brownian motion under E[-]. 

For a given p > 1, we also denote L P G (H) = {X 6 L G (H), \X\P € L G (TL)}. 
L G (Tt) is also a Banach space under the norm ||X|| := (E[|X| p ]) 1 / p . We have 
(see Appendix) 

\\X + Y\\ p <\\X\\ p + \\Y\\ p 
and, for each X e L G , Y G L q G {Q) with | + i = 1, 

\\XY\\=E[\XY\]<\\X\\ p \\X\\ q . 

With this we have < ||X|| p , if p < p' . 

We now consider the conditional expectation introduced in ([T0| . For each 
fixed t = tfc < T, the conditional expectation E[-|Ht] : LL(Ht) i— * (Tit) is 
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a continuous mapping under ||-||. Indeed, we have E[ELY|7i t ]] = ELY], Y e 
L^ p (Wt) and, since P t is subadditive, 

E[X\H t ] - E[Y\H t ] < E[Y - Y\H t ] < E[\X - Y\\H t ] 

We thus obtain 

E\E[X\Ht] - W\H t ]] < E[X - Y] 

and 

\\E[X\n t }-E[Y\n t }\\ <\\X-Y\\. 

It follows that E[-|7i t ] can be also extended as a continuous mapping Lq(Ht) >— > 
Lg(7<t). If the above T is not fixed, then we can obtain E[-|H t ] : Lq(H) 

Proposition 16 We list the properties ofE[-\Ht], t e [0, T], that hold in L® (Ht) 
and still hold for X , Y G L^(H T ): 

(i) E[X\Ht] = X, for X G L&Ht), t < T. 

(ii) IfX>Y, then E[X\H t ] > E[Y\H t }- 
(Hi) E[X\H t ] - E[Y\H t ] < E[X - Y\H t }. 

(iv) E[E[X\H t ]\H s ] =E[X\H tAs ], E\E[X\H t ]] = ELY]. 

(v) E[X + v \Ht] = E[X\H t ] + r? ; r? G L&Ht) 

(vi) E[nX\n t ] = 7]+E[X\n t } + r}-E[-X\H t ], for bounded n G L^(Ht). 

(vii) We have the following independence: 

E[X\H t ] = ELY], MX e Lg(Wt). vt > 0. 

where Lq{TL\.) is the extension, under ||-||, o/ {Tij,) which consists of random 
variables of the form , P* 2 , • • • ,B*J, G Zip(R m ), ti,-- - , t m G [0,T] ; 

m = 1, 2, • • • . ffere we denote 

B\ = B t+S — B t , s>0. 

f viii ) The increments of B are identically distributed: 

E[0(P* i , P t 4 2 , • • • , B* m )] - E[^(P tl , B t2 , • • • , B tm )]. 

The meaning of the independence in (vii) is similar to the classical one: 

Definition 17 ^4n K™ valued random variable Y G (Lq(7Y))™ is said to 6e 
independent of Ht for some given t if for each <f> G we /lave 

E^yjiWt] = E[0(r)]. 

It is seen that the above property (vii) also holds for the situation X G 
LqCH*) where Lq(H*) is the completion of the sub-lattice L>t>oLq(Ht) under 

INI- 

From the above results we have 



11 



Proposition 18 For each fixed t > 0, (-B*) s >o is a G-Brownian motion in 
Lq(H*) under the same G- expectation E[-\. 

Remark 19 We can also prove, using Lemma\^ that the time scaling of B, 
i.e., B = (vAB{/x)t>o also consists a G-Brownian motion. 

The following property is very useful 

Proposition 20 Let X, Y £ L%(H) be such that E[Y\H t ] = -E[-Y\H t ], for 
some t £ [0, T] . Then we have 

E[X + Y\H t ] = E[X\H t ] + E[Y\H t }. 

In particular, ifE[Y\H t ] = E[-Y\H t ] = 0, then E[X + Y\H t ] = E[X\H t \. 

Proof. It is simply because we have E[X + Y\H t ] < E[X\H t ] + E[Y\H t ] and 

E[X + Y\H t ] > E[X\H t ] - E[-Y\Ht\ = E[X\H t ] + E[Y\H t ]. 

■ 

Example 21 From the last relation of Provosition XTT\ -(ii), we have 

E[(AB t , B t )\ = <J A t = 2G(A)t, VA e S d . 

More general, for each s < t and n = (rf^)f , =1 G Lq(TC s ; Sd), 

E[{r)B s t ,B s t )\H s ] = a v t = 2G(n)t, s,t>0. (11) 

Definition 22 We will denote, in the rest of this paper, 

Bf=(sL,B t ), for each a= (a ir ■ ■ ,a d ) T £ R d (12) 

From Lemma\^ and Remark^ 

E[0( J Bf)]=P t G (0((a,.)))=P^(0) 

where P Gb - is the (1 -dimensional) G a -normal distribution. Thus, according to 
Definition \13\ for d-dimensional G-Brownian motion, B a forms a 1-dimensional 
G & -Brownian motion for which the G ^-expectation coincides with E[-]. 

Example 23 For each < s — t, we have 

E[i>(B t - B S )\H S ] = E[^(B t - B s )} 

If cj) is a real convex function on R and at least not growing too fast, then 

E[X<P(B*-Bf)\Ht\ 

= X+E[4>{B% - Bf)\H t ] + X-E[-4>{B% - Bf)\H t ] 

' 'x)exp( r )dx 



^2n{T-t)cr aBL T J-oo " v 2(T-t)a 



aa 1 



X~ f°° x 2 

'(x)exp(— — — - -)dx. 



y/MT-t)\a—**\J-<*> 2(T-*)|o-_ BaT |' 
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In particular, for n = 1, 2, • • • , 

n\B?-Bfr\H s ] = n\B?_ s n 



i r 00 x 2 

\x\ n exp(-— r )dx. 



^2n(t - s)cr aa T J-oo 2 (i - s)cr aa T ' 

5uf we We E[-|-Bf - B a \ n \H s ] = E[-\B?_ s \ n ] which is when a_ aa r = and 

-1 Z" 00 a; 2 

— ; / |a;|™exp(— — - :)dx, if <T- aa T < 0. 

^Mt^s)\o~~r\J-J ' Py 2{t-s)\a_ aaT \> 

Exactly as in classical cases, we have E[Bf — Bf\7i s ] — and 

E[(Bf - B») 2 \H S ] = a aa r(t - a), E[(B? - Bf) 4 \H s ] = 3a 2 aT (t - S ) 2 , 
E[(B? ~ Bff\H s ] = 15*^ (t - sf, E[(B? - Bff\H s ] = 105^ (t - s) 4 , 



a^i.o^tC*-*)^] 673 



E[|fl t a -B s a | 5 |W s ]=8 

V 71 " 

Example 24 For eac/i ra = 1,2, • • • , < i < T and X e L G (Ht), we have 
E[X(B^ - B$)\Ht] = X + E[(B^ - B?)\H t ] + X~E[-(B^ - B?)\H t ] = 0. 
This with Proposition \2(A yields 

E[Y + X(B* - B?)\H t ] = E[F|Wf], Y G L G (H). 

We also have, 

E[X(B* - B*) 2 \H t ] = X + E[{B* - Bff\H t ] + X-E[-{B* - B?) 2 \H t ] 
= [X + o- aSL T+X-a_ aSL T](T-t). 

Remark 25 It is clear that we can define an expectation E[-} on L® (H) in the 
same way as in Definition \13\ with the standard normal distribution P® (•) in the 
place of P G (-). If Id S r, then it follows from Q) that PP(') is dominated by 
Pf(-) in the sense 

p?tt)-p?W<i¥(<i>-il>). 

Then E[-} can be continuously extended to L G (Ti.). E[-} is a linear expectation 
under which (B t )t>o behaves as a Brownian motion. We have 

- E[-X] < E°[X] < E[X], ~E[~X\H t ] < E°[X\H t ] < E[X\H t }. (13) 

More generally, if V C T, since the corresponding P' = P G r' is dominated by 
P G — P Gt , thus the corresponding expectation E' is well-defined in L G (Tt) and 
E' is dominated by E: 

E'[X] -E'[Y] <E[X-Y], X,YeL G (H). 
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Such kind of extension through the above type of domination relations was dis- 
cussed in details in '441- With this domination we then can introduce a large 
kind of time consistent linear or nonlinear expectations and the correspond- 
ing conditional expectations, not necessarily to be positive homogeneous and/or 
subadditive, as continuous functionals in Lq(1~1). See Example \40\ for a further 
discussion. 

Example 26 Since 

E[2Bf(Bf - Bf)\H a ] = E[-2Bf(Bf - Bf)\H a ] = 0, 

we have, 

E[(Bf) 2 - (Bf) 2 \H s ] = E[(Sf - Bf + Bf) 2 - (Bf) 2 \H s ] 

= E[(B£ - Bf) 2 + 2{Bf - Bf)Bf\H s ] 

= ff aa T (* ~ S) 

and 

E\((Bf) 2 (Bf) 2 ) 2 \H s ] = E[{(B? - Bf + Bf) 2 - (Bf) 2 } 2 \H s ] 



= E[{(B? - Bf) 2 + 2(Bf - Bf)Bf} 2 \H s ] 

= E[(Bf - Bf) 4 + A(B? - BffBf + A{B? - Bf) 2 (Bf) 2 \H s ] 

< E[(Bf - Bf) 4 ] + AE[\Bf Bf\ 3 ]\Bf\ + Aa aaT (t s)(Bf) 2 

= ^LAt s) 2 + %^[o aaT (t - s)f/ 2 \Bf\ + Aa aSL r(t - s)(Bf) 2 



5 Ito's integral of G— Brownian motion 

5.1 Bochner's integral 

Definition 27 For T G R+, a partition ttt of [0,T] is a finite ordered subset 
7r = {ti, • • • , ijy} such that = to < ti < ■ ■ ■ < tpf = T . 

h(tt t ) = max{|tj+i - ti\, i = 0, 1, • • • ,N - 1}. 

We use 7T^ = {tg < < ■ ■ ■ < t^} to denote a sequence of partitions of [0,T] 
such that limTv^oo /j(tt^ ) = 0. 

Let p > 1 be fixed. We consider the following type of simple processes: for 
a given partition {to, • • • , iAr} = ttt of [0, T], we set 

JV-l 

Vt(u) = &( w ) I [**,*«.+o( < ) 

fc=0 

where G L^Ti^), k = 0, 1, 2, • • • , N — 1 are given. The collection of these 
type of processes is denoted by Mq°(0,T). 



14 



Definition 28 For an r, G M*'°(0,T) with n t = ££ =0 aMI [tjbl t fe+1 )(t) ; the 
related Bochner integral is 



T N-l 

/ rj t (uj)dt = V" £ k (uj)(t k+1 - t k ). 
Jo k=0 



Remark 29 We set, for each -q G Afg°(0,T), 



1 /- T 1 w_1 

M := 7^ / E[7? t ]dt = — E&M(*k+i - **) 
J Jo 1 k=0 



It is easy to check that : M G (0, T) i — ► K forms a nonlinear expectation 
satisfying (a)-(e) of Definition]]^ We then can introduce a nature norm 



\r ] \\ l T = %r[H] = ^ f nr)t\]dt 



Under this norm M G (0, T) can extended to M G (0, T) which is a Banach space. 

Definition 30 For each p > 1, we denote by M G (0,T) the completion of 
Mq (0, T) under the norm 



(r I llkl p ll*) 1/p = [yE E [I^HI p ](^+i-^ 



We observe that, 

E[| / n t {u)dt\] < V ||e fe H|| (tjk+i - **) = / 
- 70 fc=o 70 



E[|»jt|]<ft. (14) 



We then have 



Proposition 31 The linear mapping J Q rj t (uj)dt : M G °(Q,T) t— > L g (Ht) is 
continuous and thus can be continuously extended to M G (0,T) i— > L g (TIt)- We 
still denote this extended mapping by rj t (u>)dt, rj G Mq(0,T). 

Since Mg(0,T) C M^(0,T), for p > 1. Thus this definition holds for rj G 
Mg(0,T). 

5.2 Ito's integral of G— Brownian motion 

We still use Bf := (a, B t ) as in (fT2j) . 



Definition 32 for eac/i 77 G M G °(0, T) with the form r]t(ui) — J2k=o £k(u)I[ tk j k+ 
we define 

T N-l 

I{rf) = / n(s)dBf := £ - S?J. 

- 70 fc=0 
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Lemma 33 The linear mapping I : Mq (0,T) i — ► Lq(Ht) is continuous and 
thus can be continuously extended to I : M G (0, T) i — ► L g (TIt). In fact we 
have, for each r\ € Mq (0, T), 

E[/%( S )dB a ]=0, (15) 
Jo 

E[( / T](s)dBf) 2 ) < a^r f E[ V 2 (s)}d S . (16) 
Jo Jo 

Definition 34 We define, for a fixed r\ G Mg(0, T) the stochastic calculus 

pT 

V (s)dBf ;= I{rf). 

It is clear that {W\), {W\) still hold for r) G M G (0,T). 
Proof of Lemma 1331 From Example [Ml for each k, 
Efe(B» +1 -B»)|W t J=G. 

We have 

E[ [ T v(s)dBf] = E[ V (s)dB* + £, N ^{B? N - B^_J] 

JO JQ 

pt]S[ — -L 

= E[ / TjWdB' + E^-x^-^VJIWt^.J] 
Jo 

rtN—i 

= E[ v(s)dBf]. 
Jo 

We then can repeat this procedure to obtain (fT5j) . We now prove ([I 

„ 9 

ft TV- 



E[^ »7(*)dB s »j ] = E[ r?( S )dB a + Ov-i^ - Bf N i ) 

= E[ QT"" 1 7,(^5*) 

+ E[2 (7 »?(«)dB^ 6v-i(B t a „ - Bf^J + &_ x (iJ? w - Bf^flHt^}] 

ftjV-l \ 2 

J r}(s)dB*J +^_ 1 a &aT (t N -t N - 1 )\. 



Thus E[(/ * N r,(s)dBf) 2 } < E[ (j^ ilW^j'+C-J^Cff-tw-O. We 
then repeat this procedure to deduce 

E[( / j,( a )dB s ) 2 ] < o-^t V n^k) 2 ](t k+ i - t k ) = [ E[(n(t)) 2 }dt. 
Jo k=Q Jo 
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We list some main property of the Ito's integral of G-Brownian motion. Wc 
denote for some < s < t < T, 




We have 

Proposition 35 Let rj,6 e M^(0,T) and let < s < r < t < T. Then m 
Lq(Ht) we have 

(i) jlvudBl = /; Vu dB a + f* Vu dB*. ^ 

(ii) J s (ar] u + 9 u )dB a = a J g r/ u dB a + J s 6 u dB a , if a is bounded and in Lq(H. s ), 
(in) ELY + X T VudB*\H s ] = E[X\H S ], VX e L^H) 

(w) nUr^dBlf\H s ] < a aaT tfn\ Vu \2\H s }du 

5.3 Quadratic variation process of G— Brownian motion 

We now consider the quadratic variation of G-Brownian motion. It concen- 
trically reflects the characteristic of the 'uncertainty' part of the G-Brownian 
motion B. This makes a major difference from the classical Brownian motion. 
Let 71"^, N = 1, 2, • • • , be a sequence of partitions of [0, t]. We consider 



AT-l 



(B?r = j2m» +1 ) 2 -(B?») 2 } 

k=Q 

N-1 N-1 

= > 2B a N (B a N — B a N ) + > (B.N — B.n) 
^ f k v *fc+i f k ' ^ v *fc+i *fc ' 



fc=0 k=0 
r N\ _ m „ T UN +N\ 



As n(iTt ) = max o<k<N-i(tk+i — ) — ^ 0, the first term of the right side tends 
to J Q BfdBf. The second term must converge. We denote its limit by (B a ) 



i.e., 



v 

(B a ) t = lim £ (B% B^f = (Bff 2 f BfdBf. (17) 

By the above construction, {B a ) tl t > 0, is an increasing process with {B a ) = 0. 
We call it the quadratic variation process of the G-Brownian motion B a . 
Clearly (B a ) is an increasing process. It is also clear that, for each < s < t 
and for each smooth real function ip such that ip{(B a ) t _ s ) e ig(Ht- s ) we have 
Em(B a ) t _ s )] = E[i;((B a ) t - (B a ) s )]. We also have 

(B a ) t = (B- a ) t = (-B a ) t . 

It is important to keep in mind that (B a ) t is not a deterministic process unless 
the case cr aa T = — er_ aa T and thus B a becomes a classical Brownian motion. In 
fact we have 
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Lemma 36 For each < s < t < oo 

E[(B a ) t -(5 a } s |H s ]=^(i-«), (18) 
E[-((5 a ) t - (B*) S )\H S ] = «r_ aa .(t - s). (19) 

Proof. By the definition of (B a ) and Proposition [351- fiii) . then Example [26l 

E[(£ a ) f - (B a ) s = E[(B t a ) 2 - (B s a ) 2 - 2y B a <LB a |H s ] 
= n(B?f-(B*) 2 \H s }=* aa T(t-s). 

We then have (fT8"|). (|19p can be proved analogously by using the equality 
E[-((S t a ) 2 -(i? a ) 2 )|K]-a_ aa T(t- S ). . 

An interesting new phenomenon of our G-Brownian motion is that its quadratic 
process (B) also has independent increments. In fact, we have 

Lemma 37 An increment of (B a ) is the quadratic variation of the correspond- 
ing increment of B a , i.e., for each fixed s > 0, 

(B a ) t+S (B a ) s = ((B s T)t 
where B s t = B t+S - B s , t>0 and (B s )f = (a, B l s ). 

Proof. 

(B a ) t+S (B*) s = (Bf +S f — 2 J* +S B a dB a - ((Bff - 2 J* B^dB*) 
= (B? +s - B a f -2^ + \bI - Bf)dB a 



= (Bf +S - Bf) 2 - 2 J (B a +U - B a )d{B a u - Bf) 
= ((BT) t - 



Lemma 38 We have 

E[(B a f t ] = E[((B a ) t+s - (B a ) s y\H s ] = a^-rt 2 , s,t>0. (20) 
Proof. We set <j>(t) := E[(B a ) 2 ]. 



<f>(t) = E[{(B t a ) 2 - 2 / BtdBlY] 
Jo 

<2E[(B a ) 4 ]+8E[(/ BldBlf] 
Jo 

< 6al aT t 2 + 8a aa T f E[(5 a ) 2 ]d M 
Jo 
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This also implies E[((£ a ) t - (B a ) s ) 2 } = tf>(t - s) < Worrit - s) 2 . For each 

«e[o,t), 

^)=E[((B a ) s + (B a ) t -( J B a ) s ) 2 ] 

< E[((i? a ) s ) 2 ] + E[((i? a ) t - (B a ) s ) 2 ] + 2E[«i? a ) t - <5 a ) s ) (B a )J 
= <f>(s) + <f>(t -s) + 2E[E[((i? a ) t - (B*) a )\H a ] (B a ) J 
= ( />(s) + ^(t-s)+2a 2 eLBlT s(t-s). 

We set <5at = t/N, = kt/N = fc<5jv for a positive integer N. By the above 
inequalities 

0(4) < + 0(M + 2^ aT ^_ 1 5 JV 



We then have 

AT-l w-i 
fe=0 fc=0 

Let AT -> oo we have 0(t) < 2er 2 aT / Q * sds = cr^t 2 . Thus E[(S a > 2 ] < cr 2 aT i 2 . 
This with E[(5 a ) 2 ] > E°[(B a ) 2 } = cr^i 2 implies I®. In the last step, the 
classical normal distribution Pf, or 7V(0, 707J ), 70 € T, is chosen such that 

MtotJ aa T ] = cr 2 aT = supir[77 T aa T ]. 

7er 



Similarly we have 

E[((S a ) t -(B a ) s ) 3 |W s ]=a^(t- S ) 3 , (21) 
E[((B*) t -(B*) s r\H s }=*i aT (t-s)*. 

Proposition 39 Let < s < t, £ e L^(H a ), X e L^U). Then 

E[X + e((Sf) 2 - (Bf) 2 )] = E[X + - B a ) 2 ] 

= ELY + £(( J B a ) ( -( J B a )J]. 

Proof. By (fT7|) and applying Proposition [20l we have 

E[X + £((i? a ) 2 - (5 a ) 2 )] = E[X + a(B a ) t - (B a ) s + lf s B u dB u)} 

= E{X + a(B a ) t -(B*) s )}. 
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We also have 

ELY + £((Bf) 2 - (Bf) 2 )} = E[X + - Bf) 2 + 2(fl? - Bf)Bf}} 

= E[X + C(Bt-Bf) 2 }. 



Example 40 We assume that in a financial market a stock price (St*)t>o is 
observed. Let Bt — log(St), t > 0, be a 1- dimensional G-Brownian motion 
(d = 1) with T = [(T*,cr*] 7 with fixed er* G [0, \) and a* e [l,oo). Two traders 
a and b in a same bank are using their own statistics to price a contingent 
claim X = (B) T with maturity T. Suppose, for example, under the probability 
measure P a of a, B is a (classical) Brownian motion whereas under Pf, ofb, ^B 
is a Brownian motion, here P a ( resp. V^) is a classical probability measure with 
its linear expectation E a (resp. E b ) generated by the heat equation dtu = ^d 2 x u 
(resp. dtu — \d 2 x u). Since E a and E b are both dominated by E in the sense of 
(3), they can be both well-defined as a linear bounded functional in Lq(H). This 
framework cannot be provided by just using a classical probability space because 
it is known that (B} T = T, ¥ a -a.s., and (B) T — -j, P b -a.s. Thus there is no 
probability measure on with respect to which P a and P b are both absolutely 
continuous. Practically this sublinear expectation E provides a realistic tool of 
dynamic risk measure for a risk supervisor of the traders a and b: given a risk 
position X G L^Ht) we always have E[-X\H t ] > E a [-X\H t ] VE b [-X\H t } for 
the loss —X of this position. The meaning is that the supervisor uses a more 
sensitive risk measure. Clearly no linear expectation can play this role. The 
subset r represents the uncertainty of the volatility model of a risk regulator. The 
lager the subset T, the bigger the uncertainty, thus the stronger the corresponding 
E. 

It is worth to consider to create a hierarchic and dynamic risk control system 
for a bank, or a banking system, in which the Chief Risk Officer (CRO) uses 
E = E G for his risk measure and the Risk Officer the ith division of the bank 
uses E l = E Gi for his one, where 

G{A) = \ sup tr[jj T A], Gi{A) = \ sup tr^A], T.cT, i = 1, • • • , I. 

Thus E* is dominated by E for each i. For a large banking system we can even 
consider to create E 4 - 7 = E Gi ^ for its (i,j)th sub-division. The reasoning is: in 
general, a risk regulator's statistics and knowledge of a specific risk position X 
are less than a trader who is concretely involved in the business of the product 
X. 

To define the integration of a process r/ <E Mq(0, T) with respect to d (B a ), 
we first define a mapping: 

Qo,T(l) = f v{s)d{B*) s := £ U(B a ) tk+ -(B*) tk ) : m£°(0,T) ~ L\H T ). 
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Lemma 41 For each 77 <E M^ a (0,T), 

n\Qo, T (v)\] < ^ / E[\v.\]ds, (22) 
Jo 

Thus Q a ,T ■ M^°(0,T) i-> L^Wt) 

is a continuous linear mapping. Conse- 
quently, Qo,t can 6e uniquely extended to Mq(0,T). We still denote this map- 
ping by 

f V (s)d(B a ) s = Q , T (v), veMh(0,T). 
Jo 

We still have 

Efl [ V (s)d{B*) s \] <a aaL T [ T E[\ Vs \}ds, V77 S Mq(0, T). (23) 
J Jo 

Proof. By applying Lemma |3T>I (|22p can be checked as follows: 

N-l N-l 

fc=0 k=0 

N-l 

= E E [l^l] CT aa-(ifc+l-<fc) 
fc=0 

= cr aa T / E[|ry s |]ds. 

JO 

■ 

We have the following isometry 
Proposition 42 Let r) € M|(0,T) 

E[( / T «( S )di? s a ) 2 ] = E[ f r?{s)d (B a )J (24) 
Jo Jo 

Proof. We first consider r\ € M G ' (0, T) with the form 

JV-l 

m(p) = E feMWu-oW 



fc=0 



and thus J^ri(s)dBf := J2k=o &(B? h+1 ~ Bf k ). By Proposition [20] we have 

E[X + 2t k (B? k+i - B? k MB? i+i - B-)] - E[X], for X € L&H), I ? k. 
Thus 

-T /JV-l \ 2 JV-l 

E[( / V {s)dB*f] = E[ E &(Sf fc+1 - BfJ ] = E[ E ^(5f fc+1 - B fJ 2 ] 
Jo \k=o / fe=0 
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This with Proposition it follows that 



rp jy 1 rp 

E[( „( a )dBJ) 2 ]=E[££«i?»> tfc+1 -(S a ) t J]=E[/ 7? 2 ( S )d(5 a )J. 
- 70 fe=0 - 70 

Thus (HI holds for 77 € Af^°(0,T). We thus can continuously extend this 
equality to the case rj £ M<2,(0, T) and obtain (|24|) . ■ 



5.4 Mutual variation processes for G— Brownian motion 

Let a = (a%, ■ ■ ■ , ad) T and a = (01, • • • , dd) T be two given vectors in R d . We 
then have their quadratic variation process (B a ) and (-B a ). We then can define 
their mutual variation process by 

(B*,B*) t :=\[{B* + B*) t -(B*-B*) t ] 

= \\(B^) t -(B*-*) t ]. 

Since (B a - a ) = (B a - a ) = {-B a - a ), we see that (B a ,B a ) t = (B a ,B a ) t . In 
particular we have (B a , B a ) = (B a ). Let 7Tj , N — 1, 2, • • • , be a sequence of 
partitions of [0,t]. We observe that 

JV-l , JV-l 



fc=0 fc=0 

Thus as h{tt^) — > 0, we have 

JV-l 

lim V - i?fiv ) (B%t - Bf N ) = (B a , B a ) . 



fc=0 

We also have 



= -[{B a+a ) 2 -2 ( B a+a dB a+a - {B a - a f + 2 / B a ~ 5 dB a ~ a ] 
4 Jo Jo 



flffl? - / BfdB? - / BfdBf. 



Now for each 77 £ Mq(0,T) we can consistently define 

[ Vsd (B a , B~ a ) s = \ £ Vs d (B a +~ a ) s -\j\s (B^) t 
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Lemma 43 Let rj N E Mg O (0, T), N = 1, 2, • • • 



6e of form 



N-l 



fc=0 



i/jii/i /i(7T^ ) — > and - 
following convergence in L\ 



77 m M,i,(0, T) as N — > 00. TTien we /lave <fte 



c(Ht): 




jv ) (Bfiv — 




5.5 Ito's formula for G— Brownian motion 

We have the corresponding Ito's formula of $(X t ) for a "G-Ito process" X. For 
simplification, we only treat the case where the function $ is sufficiently regular. 
For notational simplification, we denote B l = B ei , the i-th coordinate of the 
G-Brownian motion B, under a given orthonormal base (ei, • • • , e<j) of M. d . 

Lemma 44 Let $ S G 2 (R") be bounded with bounded derivatives and {d^x^^Y/l 
are uniformly Lipschitz. Let s E [0, T] be fixed and let X — (X 1 , • • • , X n ) T be 
an n-dimensional process on [s, T] of the form 



where, for v = 1, • • • , n, i, j = 1, • • • , d, a" , rf^ and (3 v% i , are bounded elements 
ofL%(H s ) and X s = {X l s , ■ ■ ■ ,X™) T is a given W l -vector in L 2 G (H S ). Then we 



Here we use the Einstein convention, i.e., the above repeated indices fi, v, i and 
j (but not k) imply the summation. 

Proof. For each positive integer N we set 5 = (t— s)/N and take the partition 



X v t = X V B + a"(t -s) + V ^((B\Bi) t - (B\ &) a ) + ^{B\ - B{), 



have 




(25) 



7T| 



= -C*o " ,*%} = {8,8 + 5,- ,S + NS = t}. 
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We have 

N-l 

$(X t ) - $(X S ) = J2 M X t» +1 ) - ®( X t»)} 

fe=0 

= V[a x ,$(x tr )(xf N 
' — < k i k+i i k 

k=0 

+ \ \dl» x MX t * ) (X^ - X§ ) (X% +i - X% ) + 7^]] (26) 

where 

Vk = [dl, x MX t N +9 k (X tf+i -X tf ))-^„*(X tJ r )] ) ) 
with 6» fc € [0,1]. We have 

e[|% w |] - e[|[^$(^ + e k {x t » +i - X t ,)) - dl^(X t *)] 
x (X N — X N ){XY N — XY N )\] 

L k+1 L k "k+1 k 

<cE[\X t , +i -X t ,\ 3 ]<C[S 3 + 6 3 / 2 ] 

where c is the Lipschitz constant of „ =1 . In the last step we use 

Example l2lil and (|2"Tj) . Thus X^^U?^!] — > 0. The rest terms in the summation 
of the right side of ([H are £ t w + Cf with 

N-l 

t? = ^{d^iX^it^-t^ + ^dB^B^^ (B\Bi) t „) 

k=0 

+ ^{B{ N -Bi N )] + ^ x ^(X t ^^(BL , -BU){B{ N -B{ n )} 

L k + 1 L k Z k "k+1 "k v k+l L k 

and 

1 

x [^(tf +1 -ff)+^ m ((S i ,S m > tW — £? m ) „)] 

+ K(^ + i-^)+^'((^,^') tN — (B\ B' J ) N )][3 ul (B l N -B[ N ) 

l k+l l k k + l b k 

We observe that, for each u G [t^ ' ,t k+1 ) 

N-l 

E[\d x ^(x u ) - J2 ^ $ (^^) I [^,^ +1 )(")I 2 ] 

fe=0 

= E[\d x ^(X u )~d x ^(X t ,)\ 2 } 
< c 2 E[\X u -X t n\ 2 ] < C[S + S 2 }. 
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Thus J2k=o dx»®(X t N)I [t N tends to d x ^{X.) in M§(0,T). Similarly, 

JV-1 

E 9 2 x . x ^(X^)I Ktt N +i) (-) - in A4(0,T). 

Let TV — ► oo, by Lemma |4"51 as well as the definitions of the integrations of dt, 
dB t and d(B) t , the limit of £f in L^(Ji. t ) is just the right hand side of (f2"5")) . 
By the next Remark, we also have — > in Lg(7i t ). We then have proved 

Remark 45 Ira tte proof o/Cf -> m L 2 G (H t ), we use the following estimates: 
for ^ G M^°(0,T) sucft iftai = E^'^tf^.lW, < = {0 < 
to, • • • ,t N = T} with limN^ = and E^" 1 E [lf£l](*£+i ~ ^ C > 

for all N = 1, 2, • • • , we /mue E[| Efclo 1 ^ (*fe+i — ^ ) 2 ] — > ^ an ^> / or an 2Z fixed 
a, a eM d ; 

JV-1 JV-1 

e[| E ^(^ a )^ +I - ( sa )^) 2 i] ^ E Eiiefi -E[((i? a ) t; v +i - (i? a > 4 „) 2 i^ t «]] 

fc=0 k=0 

JV-1 

= E n\^\K a At^ +1 -tf) 2 -o, 



JV-1 



E^(( 5a ) t - +1 -( sa ) t f)(^+i-*f) 



fe=0 
JV-1 

< e Eiief Kd - • E[((s a ) 4 « x - <s a ),«)iH t «]] 

fc=0 
JV-1 

= E E [l^l]^(^ +1 -^) 2 -0, 

fe=0 



as iueZ/ as 



JV-1 JV-1 

eh E ^(tk+i - tk)(B% - b%)]\ < e mZWu+i - tk)n\B% , - ssr o 

' "* A; + l re * ■* fc + 1 fc 

-JV-1 



y^E E [i^iK^ + i-^) 3/2 
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and 



N-i 

1 ■ fc+l k k-\-l k 

k=0 
JV-1 

< E E[|^|]E[((B a ) tf+i - < J B a ) t «)|i?f„ +i - Bf„|] 

fc=0 
N-l 

< E E O^IM((s a U -<i?»U)Y /a E[|i& -sf„| 2 ] 1/2 

' ■ fc+1 A; fc+l fc 

fc=0 
JV-l 

= E E [i^i]^^(^ + i-^) 3/2 -o. 

A;=0 

We now can claim our G-Ito's formula. Consider 

X?=XS+ f a»ds + [ V ^d(B\B^ s + [ ffdBi 
Jo Jo Jo 

Proposition 46 Let a u , f}" 3 and rf i3 , v — 1, • • • ,n, i,j = 1, • • • , d be bounded 
processes of Mq(0, T). Then for each t > and in LQ^Yit) we have 

$(X t ) - HX S ) = j d x »$(X u )ffidBi + j d x MXu)<du (27) 

\d x ^{x u )^ + l -dl» x MXuW u l ^]d{B\BJ) u 

Proof. We first consider the case where a, r\ and (3 are step processes of the 
form 

JV-l 

Vt(u) = E &MI[t fc ,t fc+1 )(f)- 

k=Q 

From the above Lemma, it is clear that (|27p holds true. Now let 

X? N =Xg + f a^ N ds+ ( V ^ j ' N d(B i ,B j ) s + f l3 v s ^ N dBl 
Jo Jo s Jo 

where a , i] N and j3 are uniformly bounded step processes that converge to 
a, rj and f3 in Mq(0, T) as N — ► oo. From Lemma l44l 

) - $(X ) = / ^(X^/J^dS^ + / 9x„$ PC )<> N du (28) 
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Since 



E[\X^-X?\ 2 } 

< C f T {E[(a^ N a»f] +E[\r,Z> N - r#| 2 ] +E[(J3f N - ^f]}ds 
Jo 



We then can prove that, in Mq(0,T), 

d x <!>;.Y. V ;,/:"'• v d x „$(X.)ri uij 
d^ x ^(X N )/3^ N f3 vj ' N -> dl^iX.)^^ 
d x MX N )a v ' N - d x MX.)a v 
d x ^(X N )/3^ N d x »$(X.)[3^ 

We then can pass limit in both sides of ([28]) and thus prove ([27] 



6 G— martingales, G— convexity and Jensen's in- 
equality 

6.1 The notion of G— martingales 

We now give the notion of G-martingales: 

Definition 47 A process (M t )t>o is called a G -martingale (resp. G super mar ting ale, 
Gsubmartingale) if for each < s < t < oo, we have M t £ ^(Tit) and 

E[M t \H s ] = M s , (resp. < M s , > M s ). 

It is clear that, for a fixed X G Lq(H), E[X\H.t]t>o 1S a G-martingale. In 
general, how to characterize a G-martingale or a G-supermartingale is still a 
very interesting problem. But the following example gives an important char- 
acterization: 

Example 48 Let M G R, <f> = {4> l ) d l= i G M§(0,T;R d ) and r) = (v ij )i,j=i & 
Mq(0, T;Sd) be given and let 

M t = M Q + [ <PldBl+ f rtid{B\Bl) u - [ 2G(n u )du 7 t G [0, T]. 
Jo Jo Jo 

Then M is a G-martingale on [0, T]. To prove this it suffices to prove the case 
17 G M 2 ' (0,T;S d ), i.e., 

N-l 
k=0 
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We have, for s 6 [i/v-i,£iv], 

E[M t \H s ] = M s +E[ V g_ 1 ((B i ,Bl) t - - 2G(n tN _ 1 )(t - s)\H s ] 

= M s + E[ry« _ 2 {B\ - - B*)|«J - 2G(r ?tjv _ 1 )(t - s) 

= M S . 

In i/ie /as£ step, we apply the relation ill]) . We then can repeat this procedure, 
step by step backwardly, to prove the case s G [0, t/v-i]- 

Remark 49 It is worth to mention that for a G -martingale, in general, —M is 
not a G -martingale. But in the above example, when r\ = 0, then —M is still a 
G-martingale. This makes an essential difference of the dB part and the d (B) 
part of a G-martingale. 

6.2 G— convexity and Jensen's inequality for G— expectation 

A very interesting question is whether the well-known Jensen's inequality still 
holds for G-expectation. In the framework of (/-expectation, this problem was 
investigated in [3] in which a counterexample is given to indicate that, even for a 
linear function which is obviously convex, Jensen's inequality for g-expectation 
generally does not hold. Stimulated by this example, [315] proved that Jensen's 
inequality holds for any convex function under a g-expectation if and only if the 
corresponding generating function g = g(t, z) is super-homogeneous in z. Here 
we will discuss this problem in a quite different point of view. We will define a 
new notion of convexity: 

Definition 50 A C 2 -function h : M. i — > R is called G-convex if the following 
condition holds for each (y, z, A) £ M. x R x S f j: 

G(h'{y)A + h"{y)zz T ) - h'{y)G{A) > 0, (29) 

where h! and h" denote the first and the second derivatives of h. 

It is clear that in the special situation where G(D 2 u) = |Au, the corre- 
sponding G-convex function becomes a classical convex function. 

Lemma 51 The following two conditions are equivalent: 

(i) the function h is G-convex. 

(ii) The following Jensen inequality holds: for each T > 0, 

E[M#BzO)] > h(E[cf>(B T )}), (30) 

for each C 2 -function <f> such that h(cj>(BT)), 4>{Bt) € Lq(Ht)- 

Proof, (i) =Kii) By the definition u(t,x) := Pf[(j)]{x) = E[(j>(x + B t )] solves 
the nonlinear heat equation |T]). Here we only consider the case where u is a 
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C 1,2 -function. Otherwise we can use the language of viscosity solution as we 
did in the proof of Lemma [6] By simple calculation, we have 

d t h{u{t, x)) = h'{u)d t u = h'(u(t, x))G(D 2 u{t, x)), 

or 

d t h(u(t,x)) - G(D 2 h{u{t,x))) - f(t,x) = 0, h(u(0,x)) = h(<f>(x)), 

where we denote 

/(*, x) = ti(u{t, x))G{D 2 u(t, x)) - G{D 2 h(u(t, x))). 

Since h is G-convex, it follows that / < and thus h(u) is a G-subsolution. 
It follows from the maximum principle that h(Pf* (4>)(x)) < P t \h(<p))(x). In 
particular ([50)) holds. Thus we have (ii). 

(ii) =Ki): For a fixed (y,z,A) e R x R d x § d , we set <f>{x) := y + (x,z) + 
| (Ax,x). By the definition of P t G we have <9 t (P t G (»(2;))| t= o = G(D 2 (j>){x). By 
(ii) we have 

h(P t G (ct>)(x)) < P t G (h(4>))(x). 

Thus, for t > 0, 

\[h{P?{cp){x)) - h(4>(x))} < i[P ( G (M0))(.x) - h{cp{x))] 
We then let t tend to 0: 

h'{^x))G{D 2 ^x)) < G(D 2 xx h(cf>(x))). 

Since D x (f>(x) = z + Ax and D 2 x (j)(x) — A. We then set x — and obtain (|29|) . 
■ 

Proposition 52 T/ie following two conditions are equivalent: 

(i) the function h is G-convex. 

(ii) The following Jensen inequality holds: 

E[h(X)\H t ] > h(E[X\H t ]), t>0, (31) 
for each X S L l G (H) such that h(X) G Lq(H). 

Proof. The part (ii) =>(i) is already provided by the above Lemma. We 
can also apply this lemma to prove (f3"Tj) for the case X £ L? p (7i) of the form 
X = 4>(B tl , • • • , B tm — B tm _ 1 ) by using the procedure of the definition of E[-] 
and E[-|7i t ] given in Definition [TBI We then can extend this Jensen's inequality, 
under the norm ||-|| = E[| • |], to the general situation. ■ 

Remark 53 The above notion of G-convexity can be also applied to the case 
where the nonlinear heat equation (Qp has a more general form: 

— - G(u, X/u,D 2 u) = 0, u(0,x) = ip(x) (32) 
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(see Examples 4-3, 4-4 an d 4-5 in 144V - In this case a C 2 -function h : K i — ► M 
is called to be G-convex if the following condition holds for each (y,z,A) G 
K x R d x S d : 

G(y, h'(y)z, h'(y)A + h"{y)zz T ) - h'(y)G(y, z, A) > 0. 

We don't need the subadditivity and/or positive homogeneity of G(y, z, A). A 
particularly interesting situation is the case of g -expectation for a given gener- 
ating function g = g{y 1 z), (y,z ) € M X R d , in this case we have the following 
g -convexity: 

l -h"{y)\z\ 2 + g(h(y), h!(y)z) ~ h'(y)g(y, z) > 0. (33) 
We will discuss this situation elsewhere. 

Example 54 Let h be a G-convex function and let X G Lq(H) be such that 
h(X) G Lq(H), then Y t = h(E[X\H. t ]), t > 0, is a G-submartingale: for each 

8<t, 

K[Y t \H s ] = E[hmX\F t ])\T s ] > h{E[X\F s ]) = Y s . 

7 Stochastic differential equations 

We consider the following SDE driven by G-Brownian motion. 

X t =X + f b(X s )ds+ f h ij (X s )d(B i ,B^) s + f o- j (X s )dBi, t G [0,T]. 
Jo Jo Jo 

(34) 

where the initial condition Xq G K" is given and 

b,hij,o-j : E™ i ► M™ 

are given Lipschitz functions, i.e., \4>(x) — (f>(x')\ < K\x — x'\, for each x , x' G IR ra , 
4> = b, rjij and o~j. Here the horizon [0,T] can be arbitrarily large. The solution 
is a process X G Mq(0, T; K n ) satisfying the above SDE. We first introduce the 
following mapping on a fixed interval [0, T]: 

A.(F) :=: Y G M^(0,T;R") i — > M^(0,T;R n ) 

by setting A t = X t , t G [0,T], with 

A t (Y) = X + X a + f b(Y s )ds+ f h tJ (Y s )d(B\Bi) s + f a^dBl 
Jo Jo Jo 

We immediately have 
Lemma 55 For each Y,Y' G Mg(0, T; R"), we Ziawe tte following estimate: 

E[\A t (Y) - A t (Y')\ 2 } < c /*E[|y;-y;i 2 ]d s , tep.n 

Jo 

where the constant C depends only on K , T and the dimension n. 
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Proof. This is a direct consequence of the inequalities (fH|) . (fTB]) and (|23[) . ■ 

We now prove that SDE ([34]) has a unique solution. By multiplying e~ 2Ct on 
both sides of the above inequality and then integrate them on [0, T]. It follows 
that 



/ T E[|A t (y)-A t (y')| 2 ]e~ 2Ct *<C f e~ 2Ct f E[\Y S - Y' s \ 2 ]dsdt 
Jo Jo Jo 

= CJ J e- 2Ct dm[\Y s -Y^\ 2 ]ds 

= (2C)~ 1 C [ T (e- 2Cs - e- 2CT )E[|y s - Y^ds. 



(i 



We then have 



f n\A t (Y)-A t (Y')\ 2 ]e- 2Ct dt<±[ E[\Y t -Y{\ 2 ]e- 2Ct dt. 
Jo 1 Jo 

We observe that the following two norms are equivalent in Mq(0, T; M. n ) 

' T E[\Y t \ 2 ]dt ~ [ T E[\Y t \ 2 ]e- 2Ct dt. 
o Jo 

From this estimate we can obtain that A(Y) is a contract mapping. Conse- 
quently, we have 

Theorem 56 There exists a unique solution of X £ Mg(0, T; R" ) of the stochas- 
tic differential equation \31$. 



8 Appendix: Some inequalities in L p G {7i) 

For r > 0, 1 < p, q < oo, such that | + - = 1, we have 



a + 6| r <max{l,2 r - 1 }(|a| r + |6| r ), Va,6eM (35) 

|a6|<^ + ^. (36) 
p q 



Proposition 57 

E[|X + Y\ r ] < C r (E[\X\ r ] + E[\Y\ r ]), (37) 

E[\XY\] < E[\X\ p ] 1 / p ■ E[\Y\ q ] 1 / q , (38) 

E[\X + Y\ p ] 1/p < E[\X\ p ] 1/p + E[\Y\ P ] 1/P (39) 

In particular, for 1 < p < p' , we have E[\X\ P ]^ P < E[\X\ P ']^ P ' . 
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Proof. d37J) follows from (J33J). We set 



X 



Y 



E[|X|p]Vp 



eiifi?] 1 /?' 



By ([36]) we have 



E[\Z V \] 



< 




[^ + ^]<E[^]+E[^] 



p q 



Thus (J3HD follows. 



E[\X + Y\ P ] 



= E[\X + Y\ ■ IX + Y]?- 1 } 

< E[\X\ ■ \X + Yy- 1 ] + E[\Y\ ■ \X + Y^- 1 ] 

< E[|X| p ] 1/p • E[\X + Y\ { P~^ q ] l ' q 
+ E[\Y\ p ] 1/p -E[\X + Y\ (p - 1)q ] 1/q 



We observe that (p — l)q = p. Thus we have (|3"5j) . ■ 
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